Pair correlations in a finite-temperature ID Bose gas 
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We calculate the two-particle local correlation for an interacting ID Bose gas at finite temper- 
ature and classify various physical regimes. We present the exact numerical solution by using the 
Yang- Yang equations and Hellmann-Feynman theorem and develop analytical approaches. Our 
results draw prospects for identifying the regimes of coherent output of an atom laser, and of finite- 
temperature "fermionization" through the measurement of the rates of two-body inelastic processes, 
such as photo-association. 
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Recent observations of the one-dimensional (II?) 
regime in trapped Bose gases Q offer unique opportu- 
nities for extending our understanding of the physics of 
these quantum many-body systems. The reason is that 
the ID regime can be investigated theoretically by mak- 
ing use of the known exactly solvable II? models 0, 
which have been the subject of extensive studies since 
the pioneering works of Girardeau 0, Lieb and Lin- 
iger and Yang and Yang [f| (see d, 13 f° r re- 
views). With the development of experimentally viable 
Bose gases, these field theory models - once of theoretical 
relevance only - are now becoming testable in table-top 
experiments. 

The knowledge of the exact solutions to the ID mod- 
els allows us to go far beyond the mean-field Bogoliubov 
approximation. In the current stage of studies of exper- 
imentally feasible ID Bose gases, one of the most im- 
portant issues that requires such an approach is under- 
standing the correlation properties in the various regimes 
at finite temperatures. 

In this Letter we give an exact calculation of the finite- 
temperature two-particle local correlation for an inter- 
acting uniform ID Bose gas, g^ = (x) 2 ^ (x) 2 ) / n 2 , 
where ^(x) is the field operator and n = (iff* (x)i& (x)) 
is the linear (ID) density. As a result, we identify and 
classify various finite-temperature regimes of the ID Bose 
gas. Aside from this, the pair correlations are responsi- 
ble for the rates of inelastic collisional processes [9|, and 
are of particular importance for the studies of coherence 
properties of atom "lasers" produced in ID waveguides. 

At T = 0, the local two- and three-particle corre- 
lations of a uniform ID Bose gas have recently been 
calculated in Ref. 10]. Here one has two well-known 
and physically distinct regimes of quantum degeneracy. 
For weak couplings or high densities, the gas is in a co- 
herent or Gross-Pitaevskii (GP) regime with — > 1. 
In this regime, long-range order is destroyed by long- 
wavelength phase fluctuations 11] and the equilibrium 



state is a quasi-condensate characterized by suppressed 
density fluctuations. For strong couplings or low densi- 
ties, the gas reaches the strongly interacting or Tonks- 
Girardeau (TG) regime and undergoes "fermionization" 
0, 0] : the wave function strongly decreases as particles 
approach each other, and <?( 2 ' — > 0. 

At the non-zero temperatures studied here, one has 
to further extend the classification of different regimes. 
For strong enough couplings or low densities, we obtain 
the TG regime with g^ — ► not only at low tem- 
peratures, but also at high temperatures. In addition, 
and in contrast to previous T = results, we find a 
weak-coupling or high-density regime in which fluctu- 
ations are enhanced. Asymptotically, they reach the 
non-interacting Bose gas level of g^ — > 2 (rather than 
g( 2 ' — » 1), for any finite temperature T. 

The emergence of this behavior at low temperatures 
implies that one can identify three physically distinct 
regimes of quantum degeneracy: the strong-coupling 
TG regime of "fermionization" with g^ — > 0, a coher- 
ent GP regime with g^ ~ 1 at intermediate coupling 
strength, and a fully decohcrcnt quantum (DQ) regime 
with g^ ~ 2 at very weak couplings. In the GP regime, 
where the density fluctuations are suppressed and one 
has a quasi-condensate, the local correlation approaches 
the coherent level of g^ ~ 1. However, a free (non- 
interacting) Bose gas at any finite T must have g^ = 2. 
So, below a critical density- and temperature-dependent 
level of interaction strength one must have an increase in 
thermal fluctuations, until the ideal-gas level is reached 
in a continuous transition. At T = the transition is 
discontinuous and occurs at zero interaction strength, so 
that it can be viewed as a zero-temperature phase tran- 
sition. At high temperatures the GP regime vanishes 
and a decrease of coupling strength transforms the high- 
temperature TG regime directly into a classical ideal gas. 

We start by considering a uniform gas of N bosons 
interacting via repulsive delta-function potential in a ID 
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box of length L with periodic boundary condition. In the 
thermodynamic limit (N, L — > oo, while the density n = 
N/L is fixed) the system is exactly intcgrable by using 
the Bethe ansatz, both at T = and finite temperature 
0, . In second quantization, the Hamiltonian is 

H= — [ dx d x ¥ d x i> + - f tfcc (1) 
2m J 2 J 

where ^(x) is the field operator, m is the mass, and 
g > is the coupling constant. For Bose gases in 
highly elongated traps to be described by this ID model, 
the coupling g is expressed through the 3D scattering 
length a. This is done assuming that a is much smaller 
than the amplitude of transverse zero point oscillations 
Iq = ^JftJmujQ, where loq is the frequency of the trans- 
verse harmonic potential. For a positive a <C lo one has 



g = 2h 2 a/mlQ, 



(2) 



and the ID scattering length am is H 2 /mg ~ Iq/cl ^S> lo 
[l^| . The ID regime is reached if lo is much smaller 
than the thermal de Broglie wavelength of excitations 
At = (2nh 2 /mT) 1 / 2 and a characteristic length scale l c 
[l3j responsible for short-range correlations. On the same 
grounds as at T = |l(J, one finds that for fulfilling 
this requirement it is sufficient to satisfy the inequalities 
a < l < {1/n, A T }. 

The uniform ID Bose gas with a short-range repul- 
sive interaction can effectively be characterized by two 
parameters: the dimensionless coupling parameter 



7 = mg/h 2 n, 



(3) 



and the reduced temperature r = T/Td, where the 
temperature of quantum degeneracy is given by Td — 
H 2 n 2 /2m, in energy units {ks = 1)- 

For calculating the local two-body correlation g^ at 
any values of 7 and r we use the Hellmann-Feynman the- 
orem ^ij. At T = 0, it has been used for calculating the 
mean interaction energy and for expounding the is- 
sue of local pair correlations 10 . Consider the partition 
function Z = exp(-D/T) = Trexp(-i?/T) which deter- 
mines the free energy F. Here the trace is taken over 
the states of the system with a fixed number of particles 
in the canonical formalism. For the grand canonical de- 
scription, one has to replace the condition of a constant 
particle number by the condition of a constant chemical 
potential /i and add a term — uN to the Hamiltonian. For 
the derivative of F with respect to the coupling constant 
one has dF/dg = —Td(\ogZ)/dg and hence 



dF 



" /T dH/dg] = (Z,/2)(* t tf t *tf) . (4) 



Introducing the free energy per particle /(7,t) = F/N, 
the normalized two-particle correlation is: 
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(#ti$ft$i$f) _ 2m /9/(7,t) 
n 2 h 2 n 2 \ (?7 
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FIG. 1: The local correlation versus 7 at different r. 
The solid curves are exact numerical results, while the dashed 
curves represent analytic results (see text). 



We have calculated the free energy /(7, r) by numeri- 
cally solving the Yang- Yang exact integral equations for 
the excitation spectrum and the distribution function of 
"quasi- momenta" By implementing a post-selective 
algorithm that ensures that the derivative of /(7, r) is 
taken for constant n, we then calculate g( 2 ' from Eq. JSJ. 
The results of our calculations are presented in Fig. ^ 
We now give a physical description of different regimes 
determined by the values of the coupling constant 7 and 
the reduced temperature r. 

Strong coupling regime [7 > max(l,y / r)]. In the 
strong coupling TG regime the local correlation g^ re- 
duces dramatically due to the strong repulsion between 
particles, and becomes zero for 7 — > 00. In this regime 
the physics resembles that of free fermions, both be- 
low and above the quantum degeneracy temperature. 
Along the lines of Ref. to leading order in I/7 

the finite-temperature g^ can be expressed through 
derivatives of Green's function of free fermions G(x) — 
J dknp{k) exp (ifoc)/(27r), where np{k) are occupation 
numbers for free fermions. For the normalized local corre- 
lation we obtain g^ = 4 [(G"(0)) 2 - G"(0)G(0)j / 7 2 n 4 . 

In the regime of quantum degeneracy, t«1, the local 
correlation is dominated by the ground state distribution 
np{k) — 0(kp — k 2 ), where Uf = tsti is the Fermi momen- 
tum. Small finite-temperature corrections are obtained 
using the Sommcrfcld expansion: 
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T < 1 



(6) 



This low temperature result for g( 2 > has a simple phys- 
ical meaning. A characteristic distance related to the 
interaction between particles is aiu = h 2 /mg ~ 1/7", 
and fermionic correlations are present at interparticle dis- 
tances x > a\D- For smaller x the correlations practically 
do not change. Hence, the correlation g^> for large 7 is 
nothing else than the pair correlation for free fermions at 
a distance am. The latter is g^ ~ (kpam) 2 ~ l/7 2 > 
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which agrees with the result of Eq.JBJ) for r — > 0. 

In the temperature interval 1 <C t <C 7 2 the gas is 
nondegenerate, but the ID scattering length am is still 
much smaller than the thermal de Broglie wavelength A^. 
As the characteristic momentum of particles is now the 
thermal momentum kx = v2mT /H, one estimates that 



g (2) 

G(x) for the classical distribution 71f(/c)> we obtain: 



(kraio) ~ T /l ■ Calculating Green's function 
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(2) 



2t/7 2 , 1<t<7 2 , 



(7) 



which agrees with the given qualitative estimate. The 
local correlation g( 2 > is still much smaller than one and we 
thus have a regime of high-temperature "fcrmionization" . 

The results of Eqs. 10 and Q agree with the outcome 
of our numerical calculations. For r = 0.1 and r = 10, 
they are shown in Fig. ^ ( m the region of large 7) by 
dashed curves next to the solid curves found numerically 
for the same values of r. 

GP regime (r 2 < 7 < 1). In the intermediate cou- 
pling or GP regime, for sufficiently low temperatures the 
equilibrium state is a quasi-condensate: the density fluc- 
tuations are suppressed, but the phase fluctuates [ljj. As 
the phase coherence length 1& greatly exceeds the short- 
range characteristic length l c = h/ ^/mng, for finding the 
local correlations the field operator can be represented 
as a sum of the macroscopic component "to and a small 
component describing finite- momentum excitations. 
Actually, the component Wo contains the contribution of 
excitations with momenta k < fco "C 1^, whereas \&' in- 
cludes the contribution of larger k. At the same time, the 
momentum ko is chosen such that most of the particles 
are contained in Wo- This picture is along the lines of Rcf. 
, and the momentum ko drops out of the answer as the 
main contribution of W' to g^> is provided by excitations 
with k ~ l~ x 0| . The two-particle local correlation is 
then reduced to g^ = 1 + 2((^p') + (#'#')J/ti. The 
normal and anomalous averages, (W'ttf') and (#'#'), can 
be calculated by using the same Bogoliubov transforma- 
tion for W as in 3D. This gives the result that 
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where Ek = h 2 k 2 /2m, £k = y 1 E\ + 2ngE k is the Bogoli- 
ubov excitation energy, and n k are occupation numbers 
for the excitations. 

The integral term in Eq. (JHJ contains the contribution 
of both vacuum and thermal fluctuations. The former is 
determined by excitations with k ~ l~ , and at r = 
we immediately recover the zero-temperature result of 
Ref. |10| . For very low temperatures r ^ 7, thermal 
fluctuations give an additional correction, so that 

<? (2) = l-2V7/7r + 7rr 2 /(24 7 3 / 2 ), t< 7 <1. (9) 

The phase coherence length is determined by vacuum 
fluctuations of the phase and is 1$ ~ l c exp (iv/y/^y) |ll| . 



For r = 0.001 the above approximate result, shown in 
Fig. n a t intermediate values of 7, practically coincides 
with the corresponding exact numerical result. 

For temperatures r 3> 7, thermal fluctuations are more 
important than vacuum fluctuations. The main contri- 
bution to the local correlation is again provided by exci- 
tations with k ~ l~ , and we obtain, from Eq. (jSJ) 
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(2) = 1 + t/(2Vt) , 7 « t « yr 



(10) 



The phase coherence length is determined by long- 
wavelength phase fluctuations. The calculation similar 
to that for a trapped gas [13 gives Z « h 2 n/mT. The 
condition 1$ 3> l c , which is necessary for the existence 
of a quasi-condensate and for the applicability of the 
Bogoliubov approach, immediately yields the inequality 
r <C y/j. Thus, Eq. IjlOjl is valid under the condi- 
tion 7 <C t <C »/7, and the second term in the rhs of 
this equation is a small correction. One can easily see 
that this correction is just the relative mean square den- 
sity fluctuations. In the region of its validity, the result 
of Eq. (|10|) agrees well with our numerical data, and is 
shown in Fig. \T\for r = 0.001, in the region of small 7 
values. The exact results graphed in Fig. for different 
values of r show that the coherent or GP regime is not 
present for r > 0.1, in the sense that g^ as a function of 
7 does not develop a plateau around the value g^ = 1. 

Decoherent regime: At very weak couplings, 7 < 
mm(r 2 , \/t), the gas enters a decoherent regime. Both 
phase and density fluctuations are large. At small enough 
7 the local correlation is always close to the result for free 
bosons, g^ — 2. In this regime, the only consequence of 
quantum degeneracy is the quantum Bose distribution for 
occupation numbers of particles, so we can further divide 
it into a decoherent quantum (DQ) regime for r < 1 and 
a decoherent classical (DC) regime for r > 1. 

The result of Eq. (|10f) cannot be used for 7 = 0. In 
this case one has a gas of free bosons, and Wick's theorem 
leads to g^ = 2 at any r. For small r, our data in Fig. 
^show a sharp increase of g^ from almost 1 to almost 
2 when 7 < r 2 . This is a continuous transition from the 
quasi-condensate to the DQ regime 18]. Lowering the 
temperature lowers the value of 7 at which this transition 
occurs. For t = the transition takes place at 7 = 0. 
In this case it is discontinuous and can be regarded as a 
zero-temperature phase transition. 

The DQ regime can be treated asymptotically by em- 
ploying a standard perturbation theory with regard to 
the coupling constant g. Omitting the details of calcu- 
lations, which will be published elsewhere, for the local 
correlation we obtain 



/ 2 ) = 2-47/r 2 , ^/7<T<1. 



(11) 



At higher temperatures r > 1 the decoherent quantum 
regime (y^y < r < 1) transforms to the decoherent clas- 
sical regime (r > max{l,7 2 }) and g^ remains close to 
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2 (see Fig.nj. The local correlation is found in the same 
way as in the DQ regime and takes the asymptotic form 
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(2) _ 



= 2 - 7^271-/7-, T>max{l,7 2 }. 



(12) 
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The result of Eq. (|12|) remains valid for large val- 
ues of 7, provided that j 2 <C r. Here the de Broglie 
wavelength At becomes smaller than a\o and the regime 
of high-temperature "fermionization" continuously trans- 
forms into the decoherent regime of a classical gas. The 
corrections to = 2, given by Eqs. (fTTJl and lfT2*|l . 
are in agreement with the exact numerical calculations. 
In Fig. ^ the approximate analytical results are shown 
by the dashed lines next to the corresponding solid lines 
found numerically, for r = 0.1 and r = 1000. 

In conclusion, we have calculated the two-particle local 
correlation q( 2 ) for a uniform ID Bose gas at finite tem- 
peratures |l9j . Within their range of validity, the analyt- 
ical results agree with exact numerical calculations based 
on the Hellmann-Feynman theorem and the Yang- Yang 
equations. The knowledge of g&> allows one to deduce 
when the approximate GP equation often used for first- 
order phase-coherence calculations is valid. The value of 
g(2) i indicates that the correlation function factorizes, 
which is a necessary condition for using the GP approach. 
The prediction of coherent behavior only for certain tem- 
peratures and interaction strengths, r 2 < 7 < 1, may 
be an important criterion for atom lasers where spatial 
coherence is a necessary ingredient in obtaining interfer- 
ence and high-resolution interferometry. Our results are 
also promising for identifying the regime of "fermioniza- 
tion" in finite-temperature ID Bose gases through the 
measurement of inelastic processes in pair interatomic 
collisions such as photo-association. Here the transitions 
occur at interatomic distances which are much smaller 
than the short-range correlation length l c and, therefore, 
the rate will be proportional to the local correlation g( 2 \ 

We also find a fully decoherent quantum regime in the 
case of very weak interactions or high densities. There 
is a continuous transition from the GP regime to the 
ideal gas limit (7 — > 0) where the gas displays large ther- 
mal (Gaussian) density fluctuations with g^ = 2 at any 
finite temperature |2fJ. As 7 is decreased towards the 
ideal gas, the GP result only holds above a certain ratio 
of interaction strength to density. Below this, the GP ap- 
proach becomes invalid and there is a dramatic increase 
in fluctuations, with g^ -> 2 as 7 -> 0. 
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